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Abstract 

After a general description of the tomographic picture for classical systems, a 
tomographic description of free classical scalar fields is proposed both in a finite 
cavity and the continuum. The tomographic description is constructed in analogy 
with the classical tomographic picture of an ensemble of harmonic oscillators. The 
tomograms of a number of relevant states such as the canonical distribution, the 
classical counterpart of quantum coherent states and a new family of so called 
Gauss-Laguerre states, are discussed. Finally the Liouville equation for field states 
is described in the tomographic picture offering an alternative description of the 
dynamics of the system that can be extended naturally to other fields. 

Keyword: Tomography, Klein-Gordon equation, Liouville equation, Gaussian states, 
Gauss-Laguerre states 
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1 Introduction 



Recently it has been shown the equivalence between the tomographic picture of quantum 
states and the various standard representations of them: Schrodinger [1], Heisenberg [2], 
Wigner [5], etc. (see for instance [S] and references therein). In this paper we try to 
extend such description to classical fields. In particular we will discuss the tomographic 
description of the real scalar Klein-Gordon field inspired by the tomographic description 
of an ensemble of harmonic oscillators. In fact classical and quantum field states are 
usually considered as classical and quantum mechanics applied to describing these states 
for systems with infinite number of degrees of freedom (the field modes). Thus a state 
of a classical free field when restricted to consider just a finite number of modes can be 
treated as a statistical ensemble of harmonic oscillators. 

This attempt will generalize the description of classical (or quantum) states in two 
directions. On one side, describing classical field states involves dealing with an infinite 
number of degrees of freedom and on the other, a covariant treatment of fields implies 
taking into the description of the state its dynamical evolution. In order to show how 
to proceed with this task we will analyze the foundations of tomography for a classical 
system with a finite number of degrees of freedom and we will extend straightforwardly 
such construction to deal with classical fields. 

The tomographic description of classical systems presented here will be directly in- 
spired by the Radon transform, so that our construction can also be considered as an 
infinite dimensional extension of the Radon transform. The transition to quantum fields 
should proceed using similar ideas in the realm of quantum mechanical systems, however 
we will leave such analysis to a subsequent paper. Classical and quantum standard descrip- 
tions of the fields are dramatically different. The states of classical modes are identified 
with probability densities and the states of quantum modes are identified with Hermi- 
tian trace-class nonnegative density operators (or density matrices). The observables in 
the tomographic picture of quantum mechanics are tomographic symbols of corresponding 
operators which are constructed by means of a specific star-product scheme [6J. The anal- 
ogous tomographic representation of classical system states by means of classical Radon 
transform [7] of the classical probability density p^u) is also available [S], [S], 

Till now the classical and quantum field states have not been considered in the tomo- 
graphic probability representations, except in early attempts [10], [H] and more recently 
when applying the quantum Radon transform to study tomographic symbols of creation 
and annihilation field operators for bosons and fermions [12] [13]. The aim of our work is 
to extend the tomographic approach to the case of quantum and classical systems with 
infinite number of degrees of freedom and to introduce for classical and quantum fields 
the tomographic probability density functionals determining their states. We will also 
find the tomographic form of the classical field Liouville equation for the tomographic 
probability density functionals. 

The paper is organized as follows. In section 2 a generalized description of the to- 
mography of classical systems inspired on the Radon transform will be described. In this 
picture the description of states of a physical system as normalized positive functionals 
on the algebra of observables of the system is paramount. Notice that such framework is 
common to both classical and quantum systems. The tomographic description of a family 
of states, similar to coherent states, for an ensemble of independent harmonic oscillators 
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will be done in section 3. Then, in section 4, and using as a guideline the results ob- 
tained for the family of oscillators before, we will discuss the tomographic picture of a 
real scalar Klein-Gordon on a finite cavity. For that we will consider the field as described 
by a countable ensemble of harmonic oscillators and a family of states similar to coherent 
states for harmonic oscillators will be analyzed. It will be shown that the tomographic 
description of such states is equivalent to the original one. The tomographic description 
of the field states for the continuum case will be discussed in section 5 following similar 
lines and finally, the Liouville field equation in tomographic form will be discussed in 
section 6. Conclusions and further perspectives of this work are given in section 7. 

2 The tomographic picture of classical physical sys- 
tems: an overview 

The states of a classical system with a finite numbers of degrees of freedom are de- 
scribed by a probabihty density pi^u) on its phase space E fl. The phase space 
carries a canonical measure, the Liouville measure //LiouviUe that in canonical coordi- 
nates {q,p),q = {qi, . . . ,qn),p = has the form d/iLiouvmc(q, p) = d"gd"p = 
dgi ■ ■ ■ dg„dpi ■ ■ ■ dp„. In the case that ^2 is a domain in M^", the classical center-of-mass 
tomogram Wcm of the state p is defined as the Radon transform of the density p and 
consists of the average of p along affine hyperplanes on phase space, i.e., 

Wem(^, tx,u)^ [ p{q,p)S{X-ti.q-u. p)d'^gd>, (1) 

where /x = {pi, . . . , pn), = {^i, ■ ■ ■ ^n) and the equation X — n q — u p = determines 
an hyperplane IT in Q. The classical center-of-mass tomogram yVcm{X, fi,v) defines a 
probability density, depending on the random variable X, on the space of hyperplanes in 
Q. The state p can be reconstructed by using the inverse Radon transform: 

piq,p)^ / W,r.{X,^^,u)e^p[i{X-n■q-u■p)]dX-|^. (2) 

where d^/zd^i/ = dpi . . . dpndui . . . dz/„. 

Previous ideas can be extended by considering with more care the role of the ob- 
servables of the system in the construction of the tomogram Wcm- The description of 
a physical system involves always the selection of its algebra of observables, call it O, 
and its corresponding states, denoted by S. The outputs of measuring a given observable 
A & O when the system is in the state p are described by a probability measure pA,p 
on the real line such that pA,p{^) is the probability that the output of A belongs to the 
subset A C M. Thus a measure theory for the physical system under consideration is a 
pairing between observables A and states p assigning to pairs of them measures pA,p- In 
this setting the expected value of the observable A in the state p is given simply by an 
integral over the real line parametrized by A: 

(A), = / XdpA,p. (3) 
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Such picture applies equally well to classical and quantum systems. Thus for closed 
quantum systems the observables are described by self-adjoint operators A on a Hilbert 
space H while states are described as density operators p acting on such Hilbert space. 
The pairing above is provided by the assignment of the measure ha,p = Tr (pii^A) where 
Ea denotes the projector- valued spectral measure associated to the Hermitian operator 
A. 

The description of a classical system whose phase space is Q can be easily established 
in these terms by considering that the algebra of observables C is a class (large enough) 
of functions on f2,and that the states of the system are normalized positive functionals 
on O, thus for instance if O contains the algebra of continuous functions on Q, states 
are probability measures on phase space. If we assume that the phase space is originally 
equipped with a measure fi, for instance the Liouville measure //Liouviiie in the case of 
mechanical systems, then we may restrict ourselves to the statistical states considered by 
Boltzmann corresponding to probability measures which are absolutely continuous with 
respect to the Liouville measure, thus determined by probability densities p(w) on Q. 
We denote such space of states by S as before. Given an observable /("ccr) on fl, the 
pairing between states and observables will be realized by assigning to the observable / 
its characteristic distribution p/(A) with respect to the probability measure p(TU')dp('C!7), 
then the probability of finding the measured value of the observable / in the interval A 
is given by: 

/ P/(A)dA, (4) 

J A 

and the expected value of / on the state p will be given by: 

(/)p = / Ap/(A)dA. (5) 

JR 

The tomographic description provided by the classical center-of-mass tomograms Wcm 
above ([T]) does not allow to cope with systems whose phase space is not of the previous 
form (for instance spin systems) and it is convenient to expand the scope of the formalism 
to make it more flexible and allow for alternative and more general pictures. Other 
tomographic pictures have been proposed for both classical and quantum systems (see 
for instance |14j, and [TB] for a description of quantum tomograms in the realm of C*- 
algebras) . 

A general tomographic picture of a classical system may be given starting with a family 
of elements in O parametrized by an index x which can be discrete or continuous. Often 
a; is a point on a finite dimensional manifold that we will denote by A4, thus x & Ai. We 
will denote the observable associated to the element x by U{x) or depending on the 
context. Given a state p of the system, the correspondence x ^ U^, allows to pull-back 
the observables Ux to Ai defining the function Fp{x) on associated to the state p(tn) 
by: 

Fp{x) = {p,U{x)) := / Ux{^)pMd^,M. (6) 

The observables must be properly chosen so that previous integral is defined. For 
instance we could have chosen = f2 as in the definition of the Radon transform above 
(P), and then consider [/^ = thus the function Fp associated to the state p{'uj) 
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will be again p{vj) itself. The original state p{'uj) could be reconstructed from Fp iff the 
family of observables U{x) separate states, that is, given p ^ p' two different states, there 
exists X & M. such that Fp{x) = {p,U{x)) 7^ {p',U{x)) = Fpi{x). Then two states are 
different if and only if the corresponding representing functions Fp are different. 

Clearly up to now, our construction does not discriminate the description of classical 
systems from quantum systems. The difference will appear only at the level of the prod- 
uct structure on the induced functions Fp, as the Wigner-Weyl-Moyal approach shows. 
Another important ingredient for the tomographic description is the Radon transform. To 
give an abstract presentation of this transform, we shall assume for the time being that 
M. is a. manifold which carries a measure, so that we can consider integrable functions on 
it and perform the corresponding integrals. 

Consider now the dual space of J-'(A^), denoted as J-'(A^)', and a second auxiliary 
space M whose points will be denoted by y G M . The space M parametrizes a certain 
subspace T'(A^) C J^{M.y . In other words, for each y & M there is an assignment 
y D{y) with D{y) G a linear functional on the space of functions on Ai. We 

obtain a map from J^(A^) to J^{M) by setting for each / G J-'(AI): 

Wf{y) = {D{y)J). (7) 

For instance suppose that M parametrizes a family of submanifolds Siy) of Q, y E Af. 
If the submanifold S{y) has the form $(q, p; Xi, . . . , Xd) = Xq, y = (Xq, Xi, . . . , Xd) 
denoting a parametrization of A/", the corresponding generalized Radon transform would 
be written as: 

>V(y)= / p(g,p)5(Xo-$(q,p;Xi,...,X„))dV>, (8) 
Jn 

which has the same form as eq. ([T]). 

When the imbedding is properly chosen, it turns out that W{y) is a fair probability 
distribution on J\f which we have constructed out of the initial state p. The aim of 
tomography is to reconstruct p out of the experimental distribution functions that we 
obtain from the measurement of the selected observables parametrized by Ai. This is the 
so called inversion formula for the Radon transform. In the case that Q = R^" and Af 
denotes as in ([T]) the space of hyperplanes, then because of the homogeneity properties of 
the Dirac distribution, we find that Wcm satisfies the condition: 

>V„n(X,/X,l.) =0. (9) 

Due to the homogeneity condition Wcm depends effectively only on 2n variables 
instead of 2n + 1 and the inversion formula works, out of the "measurements" performed 
with the family of observables {/j, ■ q + u ■ p}, (/x, u) G M^", we are able to recover p by 
means of eq. ([2]). 

As an important example of the previous discussion, we introduce another kind of 
tomographic representation of the state p{q,p), the classical symplectic tomogram defined 
as: 

» n 

- PkQk - ^kPk)dqi . . . dg^dpi . . . dp„. (10) 



d d d 

X- h At ■ 7— + ■ 7— 

oX Ofjb ou 
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Notice that we have taken M. — M^", the phase space again, and J\f — J\fi x ■ ■ ■ x Mn 
with Hk the space of hnes in R^, the phase space of each individual degree of freedom 

of the physical system under consideration. Thus, wc have obtained a joint probability 
distribution of the n random variables (Xi, . . . = X. In contrast to the center-of- 
mass case, because of the presence of n Dirac distributions, we find that the symplectic 
tomogram Wp satisfies n homogeneity conditions: 



5 d d 

Xk—— + //fc— h i>k— — 

dXk o/ik ouk 



+ 1 



(11) 



In other words, the classical symplectic tomogram yVp{X, fi, v) depends effectively only 
on 2n variables instead of 3n. In fact, one can show that the symplectic tomogram 
Wp(X, /X, u) can be transformed into the center-of-mass tomogram Wcm of the same 



state p, and vice versa. Finally, out of the ' 
observables {jimk + i^kPk}, (iJ'k, J^k) eM.^,k ■ 
means of the symplectic inversion formula: 



'measurements" performed with the family of 
= 1, . . . , n, we are again able to recover p, by 



p(q, p)= ^p(X, IX, u) exp 
where d"X = dXi . . . dX^. 



i^(Xfc - iikQk - T^kPk) 



k=l 



3 Tomograms for states of an ensemble of classical 
oscillators 

3.1 The canonical ensemble 

If we consider a family of n independent one-dimensional oscillators with frequencies 
ujk > 0, its phase space will be R^" with canonical coordinates {qk,Pk), k — 1, . . . ,n. 
The Hamiltonian of the system will be if = Ylk=i-^k, Hk{qk,pk) — \{p\ + ^^1^)- "^^^ 
dynamics of the system will be given by 

Qk = Pfe, Pk = -^Iqk, k^l,...n, (13) 

and the Liouville measure on phase space takes again the form d/iLiouviUe = dgi . . . dg„dpi . . . dp, 
Making the change of variables = Qk/^/^, Vk = \/^Pk-i the dynamics is written in the 
symmetrical form 

ik = i^kflk, Vk = -(^k^k- k^l,...n. (14) 
and the Hamiltonian becomes 

n 1 

fe=i fe=i 

The Liouville measure remains unchanged under this change of variables d//Liouviiie(Q) p) — 
d^^d^p = d"^d"'77 = d^i . . . d^„d77i . . . d77„ = d//Liouviiie(^, ^) and statistical states are 
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described by probability densities p{q,p) = pi^yf])- Liouville equation determines the 
evolution of the state: 

^P = {P,^} (16) 

where the Poisson brackets are defined by the canonical commutation relations {qk,Pi} = 
^kh {ikili} = {Pk,Pi} = 0. Notice that if ■ccr = (^,77) G i7 is a point in phase space, 
then pti'^) = p{${t),ri{t)) with {^(t),r}(t)) the solution of the equations of motion ([Hj) 
starting at zu at time t = 0. 

In particular, the Gibbs state or canonical distribution is given by Pca.n{Q, p) = /Zq 
where the normalization constant Zq is easily evaluated 

Zo = [ e-'^^('^'P)d/iLiouviiie(q,p) = / e-^''S^--'=«i+''i)dei . . . d^dr^i ...d^n (17) 

n 
k=l 

Hence for a given observable / we will have: 

(/)pcan = ^ / fit ^)e-^/'^^--'=«^^i)dei . . . d^dr^i . . . dr/„. (18) 



More detailed information will be found in 

The classical tomographic description of a state p(^, ry) will be performed by means 
of a symplectic tomogram: 

~ n 

Wp{X,fj,,u)= p{^,r])Y\S{Xk- pkCk-i^kr]k)d^i---d^ndr]i...dr]n. (19) 

We recall that here we have taken Ai = M^"", the phase space again, and N' = N'l x ■ ■ ■ x A/'^ 
with A4 the space of lines in R^, the phase space of each individual one-dimensional 
oscillator. A simple computation shows that the Gibbs state tomogram reads: 



exp 



2 



2{pl + u, 



(20) 



A interesting family of states which are the classical counterpart of quantum coherent 
states can be introduced by means of the holomorphic representation (k = '^i^k + iVk) 
of phase space, hence the phase space becomes the complex space C" with the Hermitian 
structure (C, C) = Sfc=i <^k\Ck\'^- Given a point z = {zi, . . . , Zn) G we can construct 
the distribution 

P.(C,C) = iV(^)exp 
where 

n 

N{z) = W TTWfc-^exp [-ujk\zk\^] . (22) 

k=l 



^k{ZkCk + ZkCk 



.k=l 



Pcan(C,C)|^=, (21) 
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Notice that integrating the Liouville equation for such state yields: 

pz it) = p^^t), 

with 



Zk{t) 



ZkiO). 



(23) 
(24) 



The symplectic tomographic distribution corresponding to Pz{Cy C) is a product 

n 

W,^ (X, ^, ^,z) = l[ Wl':\X,, p,, z/fc, Zk), (25) 

k=l 

(k) 

where the tomogram WpJ of a single degree of freedom is a Gaussian distribution 

Wf}{Xk,pk,Pk,Zk) = 



uJk 



2vr {pl + i^, 



exp 



^k {Xk — (/ifc, t'fc, Zk))y 



2(4 + 



of the random variable Xu, with mean value 



{Xk (/Ufc, Uk, Zk)) = Pk^ (Zk) + l^k^ (Zk) 



and variance given by 



(^k 



(26) 

(27) 
(28) 



3.2 A new class of states: Gauss— Laguerre states 

We will introduce now a family of classical states, called Gauss-Laguerre (GL) states, 
inspired on the Wigner functions of the excited states of a quantum harmonic oscilla- 
tor. These functions are only quasi-distributions on phase space, however their squares 
are related to the purity of the corresponding quantum states and are true probabihty 
distributions [T7]. Thus, the family of classical states we consider is defined as: 



(29) 



k=l 



where {m} = {mi, m2, . . . , m„} is a multi-index and 



UJk_ 

2tt 



LmAY^ek + vl) 



e 2 



(30) 



Here L^^. is the Laguerre polynomial of degree and the Gaussian exponential is the 
not normalized Gibbs state ^can|/3=i- Notice that Pcl^mki^kyVk) is a classical state on a 
bi dimensional phase space. 

The symplectic Radon transform of the state factorizes: 



WGL,{m}{X,fX,u) = llyV^GlmMk,Pk,l^k). 



(31) 



k=l 
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and we will obtain 



with 







exp 


[-f] 




1 









2 {rrik — s)\ f 2s 



22mfe y J\s J (2s)! 



2s 



(32) 



(33) 



while is the Hermite polynomial of degree 2s. The above result can be obtained as 
follows. 

First, we drop the label k and write Wm{X, /i, z/) in place of Wql m i-^k, /^fc, i^k)- Thus 



W„(X,/i,z.) = ^y"L^(|(e' + r^^))e-H«^+^^)5(X-/ie-z.r^)dedry 



(34) 



dKe 



iKX 



Lin ( ^(f + ) e-t-^^'+^^^e-^^^^^+'^^Medr/. 



Now we put -\//i2 + z/2 = r^j,,/i = r ^^cosa^y^v = r^,^ sin a^,^, and ^ = r sin ^,77 
r cos ^. Then, we recast the previous formula as 



>V™(X,/i, u) 



2tt 



di^e'^^>V^(K,/i,z/) 



(35) 



where the characteristic function of Wm, i-e. its Fourier transform Wm, is given by 



2" d^ 
2^ 



2\ n2 



e 2 e d 



(36) 



The integral over the angular variable Q^v = 6' + a^v yields the Bessel function Jq, so: 



W™(i^,/x, v) 



-— T 

e 2 Jo 



X 



X d -r\. 



The above integral can be evaluated and gives ([E], n. 7.422 2) 
W„(ir,/i,z/) = e H ; 



2 V 



1 ( Kr^^ ' 

e 2 \ yzj ^ 



, ^ m 

'—y 

22m / ^ 

s=0 



2(m-s)\ /2s 



m — s 



L 



2s 



Kr 



(37) 



(38) 



where the last line has been obtained by a well known addition formula of Laguerre 
polynomials ([18], n. 8.976 3). 

We remark that the above equation yields, by multiplication over the restored label 
/c, the characteristic function WGL,{m}(-K', M, , with K = {Ki, . . . , Kk, . . . Kn), of the 
tomogram >VGL,{m}(^, J^)- 
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Besides, as yVm{K = 0,/i, z/) = 1, we get at once the normalization property of the 
tomogram Wm{X, J^)- 

Finally, we are able to perform the last integration. The Fourier anti-transform of 
Wrn{K,fi, v) is obtained by means of the integral over y = Kr^yj^fuj ([E], n. 7.418 2): 

__ i„ (, ) e eo= [-Xy) d, ^ ^ (^^) 

So, we get the predicted expression of >Vm(-^, A^, z^)- 



4 The tomographic picture of Liouville's equation 

Finally, let us discuss the tomographic form of the evolution equation for states, Liou- 
ville equation ( TTGj) . The evolution equation in the tomographic description was recently 
obtained in [19] in relation with a relativistic wave function description of harmonic os- 
cillators. We will describe it here in the realm of our previous discussion. Notice that 
because of the symplectic reconstruction formula for a classical state (1121) we can compute: 



d_ 
dt 



pit 11, -t) 



exp 



k=l 



-W,(X,^,.,t)d"X^^, (40) 



(notice that the symplectic tomogram is computed at a given fixed time) and, on the 
other hand: 



k=l 



n 

E 

k=l 
n 

E 



dH d dH d 



_dr]k d^k d^k dr]k 



d"ud"z/ 



k=l 



(27r)2" 

d";ud"i^ 
72^ 



dH d dH d 



2n 



Wp(X,/x, iy,t) 



drjkdik d^kdrjk 

dH d dH 

^k 



Pk- 



exp 
d 



dik dXk drik dXk 



exp 



Eventually, we obtain the evolution equation for the classical tomogram Wp: 



E 

k=l 

dH 

~Wk 



dt 
dH 
dr]k 



(42) 



d 



-I -1 



dXi 



_d_ 

dfXj 



d 



1 -1 



dX, 



_d_ 

dv^ 



d 



^ - 



d 



dX, 



d[ij 



Vj 



d 



' d 



dUj 



fikT7Tr^p{X,ti,iy,t) 

OJVk 



d 



dXk 



Notice that the arguments {^j} , {rjj} of the derivatives of H, for any j, are replaced 



by the operators < — 



ax. 



d 



dXi 



, respectively. Explicitly, the operator 
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['Scl ^ defined in terms of a Fourier transform as 

-1 ^ 

/ (K) exp (iKX) dK 



d 



dX 



f{K) 
iK 



exp iiKX) dK. 



(43) 



Due to the presence of such terms, for a generic Hamiltonian H the evolution tomographic 
equation is integro-differential. In the particular instance of H given by (|T5|) . because of 
the general correspondence rule: 



d 



d 



■Wo 



•9 d 



the tomographic evolution equation takes the form of a differential equation: 



dWp{X,fx,u,t) 

dt 



k=l 



_d_ 



_d_ 



Wp{X,^i,iy,t) 



cr (^{^J'k,i^k}k > l^fc OJk-^,Vk ^-^^l ) Wp(X,/x, i/,t), 



where a is the canonical symplectic form on the linear space E 



b2n 



(44) 



(45) 



5 The tomogram of the real Klein-Gordon field in a 
cavity 

Having shown that an interesting family of states for a finite ensemble of harmonic oscil- 
lators is amenable to be described tomographically, we will discuss now the Klein-Gordon 
equation for a real scalar field v^(x) in a finite cavity onl + d Minkowski space-time. Thus 
we consider Minkowski space-time M = M^"'"'^ with metric of signature (+,—,■■■ ,— ). 
Points in space-time will be written as x = {t, x) The dynamics of the real scalar field 
(p(x) = ip{t, x) is defined by the Lagrangian density: 



with Euler-Lagrange equations: 



Considering V 



m ip we get the Klein-Gordon equation: 

iptt - + rn^Lp = 0, 



(46) 



(47) 



(48) 



with A the ci-dimensional Laplacian in M°'. As we have extensively seen, tomographic 
methods are described on phase space where conjugated variables and Poisson brackets 
are available. On this carrier space dynamical equations are described by a vector field, 
first order differential equations in time. Thus, for our Klein-Gordon equations we have 
to introduce a larger carrier space where the equations will be first order in time. The 
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transition from second order equations in time to first order differential equations in time 
may be done in many ways [20] , here we shall consider one in which the new variables will 
make the equations of motion more symmetric. We would stress that by using a specific 
splitting of spacetime into a space part and a time part we break the explicit Poincare 
invariant form but of course our description is still relativistic invariant. To proceed, we 
will consider the Cauchy hypersurface C = {0} x M.'^ and the finite cavity will be defined 
as V C C. We consider the restriction of the field to the cavity V using the same notation 
ip{x) := ip{0, x),x eV and the Klein-Gordon equation becomes the evolution equation in 
the space of fields ip{x): 

^ = -(-A + mV (49) 

Boundary conditions at the boundary of the cavity V are chosen such that the operator 
—A + is strictly-positive and self-adjoint on square integrable functions on V with 
respect to the Lebesgue measure, thus we can define the invertible positive self-adjoint 
operator B = v^^^A~+77?. We will also assume for simplicity that boundary conditions 
are chosen in such a way that the spectrum of B is nondegenerate, so that the eigenvalues 
of B will be < Ui < uj2 < ... < < • • • with eigenf unctions $fe(x), B^k{x) = 
UJk^k{x), k = 1,2, . . .. Thus equation fH9l) may be transformed into a first order evolution 
differential equation by introducing the new fields: 

^ = 51/2^ ; r^ = 5-^/V (50) 

(notice that B~^^'^ is well-defined because B is positive and invertible) and the equations 
of motion fj49l) for the field cp take the simple symmetric form: 



dt il) \-B 0) U 



(51) 



Thus the equations of motion for the Klein-Gordon field constitute an infinite dimen- 
sional extension of the dynamics of a finite number of independent oscillators ( IT^ . Using 
the Fourier expansion of the fields ^ and t] with respect to the eigenfunctions $fc of B, 
^(^) = Yl'k=i^k^k{x),r]{x) = Z)fcli^fc*fc(^)) then, the mechanical variables qk = y/uJ^ik 
and pk = Vk/ can be interpreted as position and momentum for a one-dimensional 
oscillator of frequency Uk and their evolution in time, given by eq. f ll3p . as a trajectory 
in phase space VL = M?°^. In the presence of field fluctuations we have to introduce a 
statistical interpretation to the mechanical degrees of freedom {qk,Pk) or {^k,Vk) of the 
field (p{x), thus the classical statistical description of the field whose physical meaning 
corresponds to the probability of a certain fluctuation of the field to take place, will be 
provided by a probability law p on the infinite dimensional phase space R^°°. Thus in 
the presence of field fluctuations the state of the field will induce a marginal probability 
density on each mode pkiqk,Pk) defined by, 

Pk{qk,Pk)= / p{qi,q2, ■■■,qk, ■ P2,---,Pk, ■■■)Yldqidpi. (52) 

ijtk 

Such marginal probability could be understood as a probability density for the fc-th mode 
of the field ip described by the one-dimensional oscillator with Hamiltonian Hki^^kyVk)- 
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Similar considerations could be applied to finite dimensional subspaces of modes of the 
field whose statistical and tomographic description would be made as in the previous 
section. 

The canonical or Gibbs state for the field ^p{x) is given by the probability distribution 
on the infinite dimensional phase space of the system as: 



Pcan(6,6,---;^l,^2,...) = iVexp 



k>l 



(53) 



with the normalization constant to be determined by regularizing the integral: 

oo r / 1 \ 1 ~^ 

Ik 



I 



k=l 



det ( -f3B 



(54) 



what amounts to define the determinant of the operator B by using the ^-function regu- 
larization of determinants, i.e., 



with 



det[-(3B] =exp Ci^^(O) 



k=l 



(55) 



(56) 



In other words, the canonical ensemble for the real scalar Klein-Gordon field (p{x) is 
defined as the Gaussian measure with variance C = {^/3B)~^ on ]R^°°. Notice that 



1 °° 1 1 



(57) 



fe=i 



with denoting the potential U [99] of the Klein-Gordon field in the Hamiltonian 

picture. Observe that U [ip] can also be written as: 



U 



-\\B^\\^ = -{^,B^^) = - / ^{x){-A + m^)^{x)d''x. 



(58) 



Then the canonical ensemble for the Klein-Gordon field at finite temperature will be 
written in the usual form: 



(59) 



with T?Lp = Hfcli ^Qkdpk- Moreover, if F [ip] denotes an observable on the field (like the 
energy, momentum, etc.), then the expected value of F on the canonical distribution at 
temperature P will be given by: 



Vip 
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(60) 



The tomographic description of the states of the Klein-Gordon field will be performed as 
in the case of an ensemble of harmonic oscillators in section H] by choosing the space Ai 
the phase space R^°° itself and A/" = IlfcLi -^k with A/fc the space of straight lines on the 
phase space of the one-dimensional oscillator {C,k,f]k)- Then, as in f|T9|l . we will define: 



/oo 
Pcan [$,, r)] Y[ S{Xk - fik^k - J^kVk)d^kdrik (61) 
k=i 

e-m^'^i^S [X{x) - fi{x)^ix) - u{x)7]{x)] ViVri 



Here the Dirac functional distribution must be understood as an infinite continuous prod- 
uct: 

5 [X{x) - n{x)^{x) - u{x)r]{x)] = Y\_^ i^k - fik^k - J^kVk) (62) 



exp 



i / K{x) {X{x) — fi{x)^{x) — iy{x)rj{x)) d'^x 



where X{x), fi{x) and z/(x) are fields whose expansion on the modes cuk of the field (p{x) 
are given respectively by: 

oo oo oo 

X{x) = ^Xk<l>k{x); /i(x) = ^/ifc$fc(x); z/(x) = ^ i/fc$fc(a;). (63) 

k=l k=l k=l 

Notice that the time dependence of the various fields is encoded in the coefficients of 
the corresponding expansions. Taking advantage again of the scaling property of the 
delta function we may use the natural parametrization of optical tomograms defined by 
the reparametrization jlk = A^fc/A/7'|~+^ = cosOk, fjk = i^k / \/ fJ-l + ^1 ~ sin6'fc, Xk = 



Xk/ vA'|+^ and after standard computations we get: 

W;fjX, 6) = Ne-^^-^^l = iVe-/v Wd'^x ^ n^-Wx^ _ (54) 

with 6{x) = tan~^ [v{^)/^{^)] ^-nd the normalization constant N defined by choosing a 
proper regularization of the trace of the operator B. 



6 Tomographic picture of continuous modes 

If we consider the scalar field in an infinite volume cavity or in the full Minkowski space- 
time for instance, many or all of the modes of the system will become continuous. For 
simplicity we will assume that we are discussing the field in the d + 1 Minkowski space- 
time W^^^ and the continuous modes of the fields ip{x), C,{x), rj{x) are described by the 
wave vector k, say, 

= I iCke-''- + ^.ke'"-) d\ (65) 

etc. Now a state of the field ip{x) will be represented by a probability measure p [^,1]], 
again nonnegative and normalized. An example of such state will be given by the canonical 
ensemble, this is the Gaussian measure whose covariance is the operator B as in f l59|) : 

rf/ican = e-'^^t^lPy) = e-^^'^^^^^V^Vr] (66) 
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with the normahzation constant absorbed in the definition of the measure. 

We will consider as analogue of Gibbs states, states that are absolutely continuous 
with respect to the canonical state, i.e., states of the form: 

P M = / [^{x),V{x)] /^can (67) 

with 

f[ax),vix)]>0 (68) 

j f [^{x),vix)] e-^^^^^^^V^Vr] = 1. (69) 

Even though at a formal level, we may introduce as in (l6Ti) a tomographic proba- 
bility density for a state of the field of the form fl67p as a functional of three auxiliary 
tomographic fields X{x),^{x),ri{x) and apply, at the functional level, the usual Radon 
transform. The expansions ( 16 3 p will be replaced by the Fourier transform: 

X{x) = j^^ J {Xke-"'-- + X_,e"'--)d% etc. (70) 

Then, 

Wf [X{x),fi{x), u{x)] = [ f [^{x), r]{x)] 6 [X{x) - /x(x)^(x) - u{x)r]{x)] e-^^'^^'^^V^Vr]. 



(71) 

The inverse Radon transform maps the tomographic probability density given by fl7T|) 
onto the probability density functional 



K. ,] ^ J W, IX, .] e.p - M.)^ (.) - VXi.W,)V.<,) (72) 

The tomographic probability functional (ITTjl has the properties of nonnegativity and nor- 
malization, i.e. 

Wf[X{x),i2{x),u{x)]>0 (73) 

JwAx(.)M.)M.)m(.)-^. (74) 

These formulas hold true for any value of the auxiliary fields X{x),fi{x), i>{x). 

In the current case the manifold Af used to construct the generalized Radon transform 
is described by the tomographic fields X{x),i'{x), fi{x), which would be a continuum 
version of the finite-mode version of the straight lines: 

Xk - fJ'k^k -^kVk = 0- (75) 

We will end this discussion by emphasizing again the homogeneity property of the to- 
mographic description of the scalar field we just presented, homogeneity that is described 
by the condition: 



oX[x) dii[x) ovyx) 



>V/[X(a;),/i(x),z/(a;)] = (76) 
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7 The tomographic picture of evolution equation for 
classical fields 

In the previous sections we have seen that the state of the classical scalar field (p{x) can 
be described either by a probability density functional / [C,{x),rj{x)] on the field phase- 
space or by the tomographic probability density functional W/ [X{x), n{x),i'{x)]. Both 
probability density functionals are connected by the invertible functional Radon transform 
(ITTi) . (172|) and in view of this, they both contain equivalent information on the random 
field states. The dynamical evolution of states of the field ip{t,x) will be determined by 
the Klein-Gordon equation ( IFIj) 

If the Hamiltonian providing the evolution of the field is given by the sum of kinetic 
and potential energy 



n 



{^{x f + V [^{x)]) d^x 



1 



(77) 



the evolution of the probability density functional on the classical phase-space of the field 
obeys a Liouville functional differential equation: 



mx),^{y)} 



The functional Poisson brackets above are given by: 

{fB,.,I.GB.,]}^/(^KW.,fe)}^ + ^^ 

where the fields ^(x), ri{y) satisfy the relations: 

{ax),v{y)} = 5'ix-y). 
Then we obtain for the Hamiltonian above (1771) the expression: 



(78) 



my) J 



d xd y, 



dt 



,(x)^d^.- 



5^{x) 5r]{x) 



(79) 
(80) 

(81) 



which is just the n — )■ oo limit of the Liouville equation for finite number of field modes 
discussed in section HI 

In the case that V" [99] = we have the functional Liouville equation 



-^f If. ')] + / 



and the corresponding tomographic form of this equation reads 



dt 



/i(x 



dvix) 



d'^x = 0. 



(82) 



(83) 



To get this equation starting from (ITS]) we used the correspondences: 



6 



6X{x) ' 6ri{x) 
5 



^ u{x] 



SX{x) 



6jj,{x) 



SX{x) 



■,r]{x) o 



Siy(x) 



n -1 



SX{x) 



(84) 
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These relations correspond to a realization of the infinite Heisenberg-Weyl algebra 
generators (and enveloping algebra) on the field phase-space and the map of the repre- 
sentation in terms of the generator action onto the tomograms. 

The rule fl8^ provide a possibility to construct the tomographic form of the Liouville 
equation f ET]) . Using the substitution / — W/ in f ET]) . and the substitutions ( 151|) . we get 
the field evolution equation 



^Wf[x,f^,u] = J dVM 



d'^x' 



6V 



sax') 



5fi{x) 



6v{x' 
5 

5X{x)_ 



-V 



v{x) 



5X{x') 



W/[X,/i,z/].(85) 



For the case of field which is a collection of noninteracting oscillators described by the 
potential energy 

V [^] = -mV' (86) 



then 



reads 



dt 



f 



d'^x ( yu(x) 



f 



u{x) 



5W 



f 



which is the equivalent of eq. 



6i'{x) ^ ' 5^{x) 
to the continuous scalar field (p. 



(87) 



8 Conclusion and perspectives 

A proposal for the tomographic description of a family of statistical states for a classical 
real scalar Klein-Gordon field has been presented inspired by the tomographic description 
of statistical states for an ensemble of harmonic oscillators. This tomographic description 
of classical fields shares most of the tomographic properties of tomograms for classical 
states: homogeneity, positivity and normalization. Moreover the field equations, repre- 
sented as the evolution equation for field states, are reproduced in tomographic terms, 
paving the way towards a tomographic description of the quantum scalar field. Notice 
that the tomographic description presented in this work, a natural extension of Radon 
transform, breaks the Lorentz covariance of the field theory, thus the Lorentz covariance 
of the tomographic description should be restored at the end. Lorentz covariance, as well 
as gauge invariance (when interactions are introduced), should be incorporated as a nat- 
ural ingredient in the tomographic picture. The tomographic picture of other fields like 
Maxwell, Dirac, Proca, Einstein could be addressed following similar arguments. Such 
issues will be discussed in subsequent works. 
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